On the shape of vortices for a rotating Bose Einstein condensate 
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For a Bose-Einstein condensate placed in a rotating trap, we study the simplified energy of a vortex 
line derived in Q in order to determine the shape of the vortex line according to the rotational 
velocity and the elongation of the condensate. The energy reflects the competition between the 
length of the vortex which needs to be minimized taking into account the anisotropy of the trap 
and the rotation term which pushes the vortex along the z axis. We prove that if the condensate 
has the shape of a pancake, the vortex stays straight along the z axis while in the case of a cigar, 
the vortex is bent. 



PACS numbers: 03.75.Fi,02.70.-c 



I. INTRODUCTION 



Dilute Bose-Einstein condensates have recently been 
achieved in confined alkali-metal gases and the study of 
vortices is one of the key issues. One type of experiments 
consists in imposing a laser beam on the magnetic trap 
holding the atoms to create a harmonic anisotropic ro- 
tating potential ||, ||, ^, ^. Vortices are nucleated and 
the number of vortices depends on the rotational velocity. 
It has been observed experimentally that when the 
first vortex is nucleated, the contrast is not 100% which 
means that the vortex line is not straight but bending. 
Numerical computations of the Gross Pitaevskii equation 
have shown evidence of vortex bending ||^ . 

The aim of this paper is to characterize the dependence 
of the shape of the vortex line on the elongation of the 
trap and the rotational velocity. In particular, using a 
simplified energy for a vortex line derived in ||^ from 
the Gross Pitaevskii energy, we study the stability and 
instability of the straight vortex and we prove that when 
the condensate has a cigar shape the first vortex is bent, 
while when it is a pancake, the first vortex is straight and 
lies on the axis of rotation. We also show that vortices 
cannot be nucleated too close to the boundary, because 
they have a minimal length. 

In [^, we have derived a simplified expression for the 
energy of several vortex lines in a rotating trap from the 
usual Gross Pitaevskii energy describing the steady state 
of the condensate 
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We let d — {fi/mujyY/'^ be the characteristic length, 
LOx — OLOjy, LOz = P^y Wc define a small non-dimensional 
parameter e which characterizes the fact that we are in 
the Thomas Fermi regime by 



where N is the number of particles and a the scattering 
length. In the ENS experiment H, |l, e = 1.74 10"^ 
while in the MIT experiment 0, e = 3.52 10"^. We 
rescale distances by d/^ and the chemical potential /ig 
so that the new chemical potential po is given by 
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In these units, we have po = 0.42 and po = 0.46 respec- 
tively for the ENS and MIT experiments. We let 



p(r) = po-(aV+2/" + /3"z") 
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be the Thomas Fermi limit of the wave function in 
rescaled units. Then, we have obtained in Q a simplified 
expression for the energy of a vortex line 7, which is 



efiujyTT\\og e\ E['^] 



with 



p dl — fl / p dz, 
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where n is related to the experimental rotational velocity 
n by 



n 



1 



ujy {l + a'^)e\log e\ 
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The energy E[^] reflects the competition between the vor- 
tex energy due to its length (1st term) and the rotation 
term. Note that the rotation term is an oriented integral 
{dz not dl), which actually forces the vortex to be along 
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the z axis, while the other term wants to minimize the 
length. This is why, according to the geometry of the 
trap, the shape of the vortex varies. 

This energy is similar to the one obtained in 11 in the 
study of the dynamics of the vortex line. Note that the 
energy that we actually derive in |^ is slightly more in- 
volved than (^. In the regime of the experiments, it is 
reasonable to restrict to this expression (||), taking into 
account the fact that the vortex core is sufficiently small 
(it is of size e in our units) and neglecting the interaction 
of the curve with itself. We are interested only in the 
presence of the first vortex: when there are several vor- 
tices, the energy has an extra term due to the repulsion 
between the lines. 

In this scaling, the energy of the vortex free solution is 
zero. Thus, a vortex line is energetically favorable when 
ri,/3 are such that inf-yi?[7] < 0. The aim of this paper 
is to study the shape of the vortex lines 7 minimizing 
EYi\. We define the domain V = {p > 0}. This is the 
domain where the condensate lies. All the analysis will 
be made in 2?. In what follows, we assume that we are at 
a velocity such that there is a vortex line and we want 
to find conditions on 17 and the elongation /3 for the line 
to be stable and either straight or bent. 

First of all, it has been observed numerically Q that 
the vortex line lies in the plane closest to the axis of 
rotation and we can provide a rigorous justification: 

Theorem 1 //a < 1, then the energy is minimized when 
the vortex line lies in the {y, z) plane, that is the plane 
closest to the axis. 

Indeed, if we have a curve 7 parametrized as ^(t) = 
{x{t), y{t), z{t)), then we can define the new curve 7(i) — 
{Q,y{t),~z{t)) by ~z{t) = z{t) and y{t) ^ ^/a^x^Ty^. 
Then p{"i{t)) = p{l{t))- Since a < I, f < x^ + y^ , 
hence p(7)|7| — ^p{'^)z < p{'^)\j\ — ilp{'^)z. It follows 
that the energy of the new curve E['y] is less or equal 
than £'[7]. If a = 1, that is the cross section is a disc, 
then our arguments imply that the vortex line is planar, 
but of course all transversal planes are equivalent. 

From now on, we will assume that the curve lies in 
the plane {y,z), so that p, given by (H), only depends 
on y and z. Recall from the expression of E, (^), that 
for EIj] to be negative, we need p — Vlp^ to be negative 
somewhere, that is fip > 1. For fixed fi, we define the 
regions 

Vg:^{{y,z):^p{y,z)>l}, Vb:^V\Vg. 

We will refer to these sets as "the good region" Vg and 
"the bad region" Vi, respectively. In the bad region, the 
energy of a vortex per unit arc length is necessarily pos- 
itive, since p — ftp^ > 0, whereas in the good region, for 
appropriately oriented vortices it can be negative since 
p— Vlp^ < 0. One can see easily that for 7 to have a neg- 
ative energy, part of the vortex line has to lie in the good 
region, that is close to the center of the cloud. Note that 
for Vg to be non empty, we need at least ilpo > 1. In 



the region Vg, we will see that the vortex is close to the 
axis for all /3. On the other hand, in the region Vf,, the 
vortex goes to the boundary along the quickest path: if 
f3 is small, perpendicularly to the boundary, which gives 
rise to a bent vortex and if /? > 1, the vortex stays along 
the axis of rotation. 

The organization of the paper is the following: first we 
study the local stability of the straight vortex: if is 
large, then the straight vortex is a local minimizer. That 
is when gets large, the vortices tend to be straight, 
while if (3 is small then the straight vortex loses local sta- 
bility and the first vortex to be nucleated is bent. Then, 
we study the minimization of EIj] in Vg and 2?h accord- 
ing to the value of /3. We finally derive that a minimizer 
of the energy has a minimal length. 



II. STABILITY AND INSTABILITY OF THE 
STRAIGHT VORTEX 

In this section, we study the stability of the straight 
vortex. We parametrize the straight vortex as jsiz) = 

(0, z) for -Zmax < Z < Z^ax, with Zmax = y/PO / P ■ One 

can compute -B[7s] and derive that it is for fipo = 5/4. 
We have two aims: first show that for (3 small, when the 
straight vortex has energy or small negative energy, 
that is for Jlpo close to 5/4, then it is unstable. Then, 
we want to prove on the contrary that if (3 is fixed and 

is sufficiently big, the straight vortex is stable. 

We consider perturbations of the straight vortex of the 
form 75(z) = {5v{z), z + 5'^w{z)) + 0{5'^) for \z\ < z^ax- 
We require that w be chosen so that p(75(±Zinax)) = 
0, thereby respecting the condition that the vortex line 
terminate at the boundary of the cloud. 

Writing a Taylor series expansion for E, one finds that 

E[js]^E[j,] + ^-^{v,E"[j,]v) + 0{S^) 

where 

2{2np-l)v^ + pv^dz. 

Here and in the rest of this section, p = p{0, z) = po — 
P'^z'^. To get this it is necessary to integrate by parts 
and use the fact that the straight vortex solves the Euler- 
Lagrange equations for E. In particular this eliminates 
all terms involving w. No boundary terms arise from 
integration by parts because ^(75) = at the endpoints. 

We say that the straight vortex is stable if 
{v, E"['ys]v) > for aU v, and unstable if {v, E"[-fs]v) < 
for some v. 

Theorem 2 The straight vortex is stable if 

3 1 
"^°>4 + 4^- 
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The straight vortex is unstable if f3 < \ /\/3 and 



(6) 



Note that the 2 values are consistent in the sense that 
they both scale like l//3^ when (3 is small. For Q, large, 
one expects several vortices in the condensate, but the 
fact that a straight vortex is stable gives an indication 
that for SI large, each vortex should be nearly straight, 
which is consistent with the observations 1^. 

Remark 1 It is interesting to see what happens in The- 
orem 1^ when ilpo = 5/4, that is when the straight vor- 
tex has zero energy. The first inequality yields that if 
(3 > 1/^/2, then the straight vortex is stable for all fl such 
that flpo > 5/4, that is when i?[7s] < 0. If (3 > \ , we will 
see that 7^ is not just stable but in fact minimizes E. The 
second inequality implies that, if (3 < y^/lS ~ .39 then 
the straight vortex is unstable at the velocity VLpo = 5/4 
at which E[^s] — 0. As a result, for these values of (3, 
the first vortex to nucleate as increases is a bent vortex. 
Note that it has been observed in ^ that for [3 < 1/2, the 
ground state of the system corresponds to a curved line. 

All this indicates that by varying the elongation of the 
condensate, one may hope to go from a situation where 
the first vortex is bent, to a situation where it is straight. 

To prove the instability of the straight vortex, we will 
find explicit perturbations v for which {v, E"['ys]v) < 0. 
These also indicate the shape of good test functions. 

We define a perturbation v (depending on a parameter 
9, which for now we regard as fixed) by 



v{z) = 







e ui- 



i{ z < 9z^!^^ 

if Z > 0Z,„ax. 



Here v is normalized so that w(zmax) = 1- For this choice 
of V, a lengthy but straightforward calculation shows that 

3/2 

(v,E"l-,.\v) =32£!I_[(i_ 9)2,9 + 4) 
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It follows that the straight vortex is unstable if 



(1 - ef{e + 4) < — (^(1 - - + -) j (7) 

for some 6* G [0, 1). It is helpful to write 6* as 6' = 1 — rj0^ 
for some > to be determined. Then (|^) can be written 
in terms of 77, as 

^^^ l + (/3V2)-(3/2^) ^ 



This is satisfied if 

1 + (/3V2) - (3/2r,) 



flpo < 



1 

2^ 



1 



2(1 



2r] 



The extremum is achieved for rj close to 3, so we can 
take ry = 3 to find that is a sufiicient condition for 
instability. Because 6 = 1 — rj0^ > 0, this conclusion only 
holds if /3 < 1/ Vs. For larger values of /?, one can make 
different choices of to find thresholds for instability. 

To derive the sufficient condition for stability, note that 
for every z. 



2po 



2po 



1. 



Multiplying v"^ by the expression on the left and integrat- 
ing by parts, we obtain 



V dz = I p 



3v^ z , 
— + —vv 
2/30 Po 



dz. 



Since \z\/pQ < Zmax/po = ^/Py/po for \z\ < z„ 



v^dz < p 



\v\ \v'\ 



dz. 



2po 13^^ 
Now we use the inequality ab < + 6^/2, to deduce 



^max f^max 

v'^dz < 



2po 2poP 2 



dz. 



In particular, if 



3 1 

"^°>4 + 4^ 

then this implies that {v, E"['ys]v) > for all v. This 
completes the proof of Theorem |[ 

We would like to derive a more precise estimate of the 
critical velocity for which a bent vortex minimizes the 
energy £^[7]. We have seen that for £^[7] to be negative, 
we need at least ftpo > 1 so that the good region T>g is 
nonempty. Note that ^Ipo = 1 is exactly the 2D critical 
velocity at the plane z = for the existence of a vortex. 
But a bent vortex cannot be a minimizer of £[7] exactly 
at Qpo = 1, since the good region T>g has to have some 
critical size so that the vortex energy in the good region 
provides a sufficient contribution to compensate the pos- 
itive part due to the length in the bad region. On the 
other hand, for fipo = 5/4, the straight vortex has en- 
ergy. Thus, the critical velocity to obtain a bent vortex is 
1 < OcPo < 5/4. We want to obtain a sharper estimate 
by using appropriate test functions. To find good test 
functions, note that 



A'{9) = [39^+A9-{7 



2npo 



{2 + (3')] 
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and so A has a local maximum at 



25 



(1-/32) 



which lies in the interval (0, 1) for the parameter range 
that we care about. 

Note that 6^ is an increasing function of il, which is 
consistent with numerical calculations showing that for 
larger values of 17, the minimizing path stays close to the 
z axis over a longer interval. For 9 = 9^, we compute the 
energy of the path which is straight between z — —9 and 
9 and goes to the boundary along a straight line. The 
optimal end point on the boundary is at z = + /3 for (3 
small. For this special test function 7, we can compute 
E['y] to find that it is less than 

- 36/3 + 21/32 _ 4^3) _ ^ + 3^ _ ^2 ^ ^0 _ 8^^^ 
8 4 8 

Thus, for (3 small, we find an upper bound for the critical 
velocity which yields a negative energy for such a test 
function: 



(220 + 96/3) 
(203 + 76/3) ■ 



In the condition of the ENS experiment, this yields ilpo < 
1.08, that is in the original variable (see (||)), Vl/ujy < 
0.385, which is very close to the value found numerically 
0.38 §. 

As a conclusion, we have shown that there is a critical 
value of called ilc with flcPo ~ 1.08, such that a bent 
vortex has negative energy and less energy than a straight 
vortex. 



III. SHAPE OF THE VORTEX ACCORDING TO 



/3 



In this section we prove that when the condensate 
cloud has a pancake shape, then the straight vortex is 
always minimizing among vortices with negative energy. 

Recall that V = {p > 0} and we write ^{t) — 
{y{t),z{t)) to denote a generic vortex line represented 
by a continuous Lipschitz function from / = [0, 1] into V 
such that 7(0), 7(1) e dV. 

For such a curve 7, let I^,g :— {t £ I : j{t) £ Vg} and 
I-yfi = I \ 1-1, g- And let jg be the restriction of 7(-) to 
Ij^g, and similarly 7;,. 

The definition of I^^b implies that p{'-f{t)) — V,p'^{j{t)) > 
for t £ I-y,b, and as a consequence 

p{nmi{t)\ - ^p\im > m\ (pm) - ^p'm)) 

which is positive in I^^i,. Thus if 7 is such that I^^g is 
empty, then clearly £^[7] > and it is energetically fa- 
vorable not to have a vortex. This is the case in particular 
for Opo < 1 since then Vg is empty. We may thus restrict 
our attention to the case Ij_b nonempty. 



Proposition 1 Let Mg = inf{£'[7g]}, where jg is the 
restriction of ^{■) to I-y.g, then for all [3 and all VL, the 
infimum is attained by the straight vortex. 

Proposition 2 Let Mi, = inf{iJ[7;,]}, , where 7b is the 
restriction of"f{-) to I-y^b, then for all /3 > 1, the infimum 
is attained by the straight vortex. 

Note that in the bad region. Proposition 2 only holds for 
/3 > 1. If /3 < 1, the situation is somewhat more compli- 
cated: /^^ pdl is minimized by a path that joins Dg to dD 
along the y axis, whereas — J^^ p^dz is minimized by the 
straight vortex running along the z-axis. The minimizer 
of the full energy reflects the competition between these 
two terms, and hence is bent. 
We always have 

£[7] = E[-fg\ + E[-ib] > Mg + Mb 

In particular, as a corollary of the above propositions we 
deduce 



Theorem 3 For /3 > 1 

£[7] >inf(0,£;[7s]) 



(8) 



where 7s is the straight vortex along the z axis. If E['ys\ < 
0, the equality in (Q) can happen only if 7 is the straight 
vortex. 

To prove Proposition |^, first note that 
/ pdl - Vlp^dz > / p\dz\ ~ Vlp^dz > {p - D.p^)dz. 

•'Ig •'Ig Jig 

Since we have assumed that 7 does not self-intersect, we 
can identify 7 with the (oriented) boundary of an open 
set V C T). Then jg can be identified with Vg n dV — 
d{Vg nV)\ {dVg n V). Since p - flp'^ = precisely on 
dVg, this implies that 

/ {p~np^)dz = I {p-np^)dz. 
Jig JdiVgnv) 

And by Stokes' Theorem, 

{p — ^lp^)dz ~ (1 — 2^lp)pydydz. 

d{Vgnv) JvgHV 

The definition of implies that 1 — 2f7p < 0, and so this 
integral is clearly minimized if Vg n y is just the subset 
of Vg where py > 0, so that 



[ {p~np^)dz > I 



)GX)g:y<0} 



(1 — 2flp)pydydz. 

(9) 

Again using Stokes Theorem and the fact that p — ilp^ 
vanishes on dVg, we find that this is equal to 



(p(0,z)-r!p2(0,z) 



dz. 
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where (0,±z*) are the points where the z-axis intersects 
dVg. Combining these inequaUties, we find that 



f pdl-np^dz> f (p{0,z)-np^{0,z)\ dz. (10) 

7g J-^, 



It is easy to see that equahty holds in (^), and hence in 
(p^), exactly when 7 is the straight vortex, and so we 
have proved Proposition |l|. 

To prove Proposition H, fix 7 such that /^.g is 
nonempty. The beginning and end of 7 must lie in the 
bad region, and 7 intersects the good region, so I-y^b must 
consist of at least two components. Let (ai,6i) denote 
the first such component and (02, 62) denote the last, and 
write 71 and 72 to denote the corresponding portions of 
76, so that 71 is parametrized as 71 = {y,z) : (ai,6i) 
Vb, with 7i(ai) £ dV and 71(61) G dT>g. We need to 
show that 7i and 72 both have more energy than the 
corresponding parts of the straight vortex. We will con- 
sider only 71, as the argument for 72 is exactly the same. 

Define 7^ = (0, C) to be a parametrization of the part of 
the straight vortex joining (0, — Zmax) to (0, — z*), where 



~m^~-Mt?+P^m'f'^ max C(s). 

p a<s<t 

Recall that we have 71 = {y{t), z{t)). The definition is 
arranged so that 1 1-^ ({t) is nondecreasing and |7s| = C- 
To prove the proposition, it thus suffices to show that 

p(7i)|7i| - f^p2(7i)i > P(7.)l7s| - ^^p'(7s)C• 

If ({t) > ({t), this is clear, because then <^ = 0, so the 
right-hand side vanishes while the left-hand side is non- 
negative, by the defining property of the bad region Vh- 
And if C(0 = C{t), then p{'-fi{t)) = p{'-fs{t)), and so in 
this case < 1 - np{ji{t)) = 1 - np['-fs{t)) < 1. So we 
only need to show that 



I7I - ci > l-jsl - cC 

for any c G [0, 1]. We will apply it to c = Qp{'^s{t))- 
To do this, first note that 



(11) 



So 



ICI < l7l(i(|i + ^^) 



1/2 



Since (3 > 1, we conclude that \C, \ < Also, it is clear 
that |i| < |7i|. So if < a < 1, then 

I71I -cz> |7i|(l - c) > C(l - c) = |7,| - cC, 
which proves (O) , and hence Proposition 0. 



IV. MINIMAL LENGTH 

In the case /? < 1, that is when the vortex line is bent, 
we will prove that the vortex has a minimum length. 
This is related to the fact that the vortex has to go to 
the center of the cloud and spend some time in the good 
region. 

For an open set U <Z T) with Lipschitz boundary, we 
endow dU with an orientation in the standard way, so 
that Stokes' theorem holds. 

We will prove the following isoperimetric-type inequal- 
ity: 



Theorem 4 For every < /3 < 1 



p dz 



au 



< (2^^)1/2 (^^^pd/^ 



3/2 



(12) 



for every connected open subset U dV, 



Remark 2 The exponent 3/2 is the best possible. An 
inequality similar to (|7^ is valid for /3 > 1, but the proof 
needs to be modified a bit. For the straight radial vortex, 



p dz — 



idu 
and so 



16 (Po) 
15 (3 



5/2 



p^dz 



au 



pdl 



au 



-3/2 



„ 4(po)^/^ 
au 6 p 



0.52 P''\p,Y'\ 



This shows that the constant {2^^)^^"^ in ( [7^ is fairly 
close to sharp for j < P < 1 say. 

1. We use Stokes' Theorem to calculate 

/ p^dz = 2 I ppydydz < 2 ppydydz 

JdU Ju Ju- 

where U~ = {(y, z) € U : y < 0}, since ppy < for 
{y, z) £T> such that y >0. 

So the coarea formula implies that 



p'dz<2[ p}^\Vp\dydz 
au Ju- 




{iy.z)eU- : p{y,z)=s} 



|Vp| 



dl ds 



where = iiif{p{y,z) : {y,z) G U}, and p* = 
snp{p{y, z) : (y, z) e U}. Thus 




< \p* - 



{iy,z)eU : piy,z)=s} 



Py 
|Vp| 



dl 
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Thus to prove the theorem it suffices to estabhsh the 
following two claims: 



Pv 

l{(y,z)(^U : p{y.z)=s} \^ P\ 

for every s, and 



dl < 



pdl (13) 



1/2 



pdl 



dU 



dU 



1/2 



(14) 



2. We first prove (p^). Fix some s € {p„,p*) and write 
Ts to denote {{y,z) £ : p(y,z) — s}. Also, let Tg 
denote dU n {p > s}. 

First assume for simplicity that Tg is connected, so 
that it consists of the short arc of the ellipse {p — s} 
joining two points, say po = (j/c^o) and pi = (yi,zi) 
with zq < zi. We can represent Tg as the image of the 
mapping 

z ^^ {y{z), z) = {-{s - f3^z^)^/^, z), zo< z < zi. 

Differentiating the identity p{y{z),z) = s we find that 
Pyy'{z) + = 0. Thus 



d 



^{l + y'{zff''^ 



dz 

It follows that 




\Py\ ■ 



{(y.z)<^U : p{y,. 



Pv 

,z)=s} |Vp| 



dl = 



dz. 



On the other hand, the one-dimensional measure of V s is 
certainly greater than \pi — pn\ ^ zi — zq, and p > s on 
Fg, and so 



p{z,y)dl > s KTs) > s{z2-zi). 



This proves ( |13[ ) if F^ is connected. If not, one can apply 
the same argument on each connected component of F^. 

3. Next we prove (p^). Let and q* be points in dU 
such that p{q*) = p*, p{q*) — p* ■ Since we have assumed 
that U is connected, dU contains a path joining q^, to 
q* . In fact it contains two such paths. If we write V to 
denote the set of all Lipschitz paths in T) joining the level 
set {p = p*} and the level set {p = p*}, it follows that 



du 



p dl > 2 inf / p dl. 



Arguments in the proof of Proposition g show that for 
/? < 1, inf^gp pdl is attained by a path that goes in a 
straight line along the y axis. Thus 



inf p dl 



(po - y"^) dy, 



where y^ = ^po - p*, y* = VPo - P* ■ And since y* , y* < 



(po - y^) rfy > I (po - y^)2y dy 



1 /-p* 



'2^/Po Jp, 



pdp 



1 



Since 6^ — > (6 — a)^ when < a < 5, we deduce that 
(|r^) holds. This concludes the proof of the theorem. 



A short calculation starting from ( |12D shows that if 
^;[7] < then 



(15) 



We expect that even for a configuration with multiple 
vortices, adding a new vortex would require to overcome 
an energy barrier which implies that the length of the 
new vortex has a lower bound of the type (|l5|). 



CONCLUSION 



We have studied the shape of the first vortex line to 
be nucleated in a harmonic anisotropic rotating poten- 
tial, according to and the elongation of the cloud /?. 
We investigate the stability of the straight vortex and 
get that when fi is large, the straight vortex is a local 
minimum of the energy. We prove that when a vortex 
is nucleated, it is close to the axis of rotation where the 
condensate density is high, and that near the boundary, 
where the density is low, the shape of the vortex depends 
on whether the cloud has a cigar or pancake shape. This 
shape reflects the competition in the energy between the 
rotation and the inhomogeneity of the trap, which makes 
the geometry of the experiment very important. In the 
case (3 > 1 (pancake), the vortex stays straight along 
the z axis while in the case /? small (cigar), the vortex 
is bending. In the case f3 small, this allows us to define 
an energy barrier for the nucleation of vortices and to 
prove that when a vortex line is nucleated, it has to have 
a minimal length. 
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